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Energy spectrums of a nonrelativistic particle and an H-like atom in a spherical box of size R with 
general conditions of "not going out" through the box surface are explored. The lowest energy levels 
reconstruction is described from the point of view of their asymptotical behavior for large R. The 
role of von Neumann- Wigner level reflection/avoided crossing effect in this spectrum reconstruction 
is emphasized. The properties of atomic H ground state in a cell, formed by a spherical cavity with 
an outer potential shell and Neumann condition on the outward boundary, are studied in detail. 
Some of them turn out to be quite new. The relevance of such a cell to a cubic lattice of cavities, 
occupied by H, is discussed be means of first principles and assumptions of the Wigner-Seitz model. 
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1. Introduction 

Considerable amount of theoretical and experimental 
activity has been focused recently on spatially confined 
atoms and molecules [l]-[6]. So far, starting from the 
works of Michcls [7] and Sommerfeld [8], main atten- 
tion has been devoted to the properties of atoms and 
molecules, confined by an impenetrable or partially pene- 
trable potential barrier ([9], [10] and refs. therein). How- 
ever, in reality general boundary conditions of "not going 
out" don't unavoidably imply genuine trapping of a par- 
ticle by a cavity, rather they could in some special cases 
correspond to a quite different picture, where the par- 
ticle state undergoes derealization from the box with 
definite symmetry properties of the wavefunction, as in 
the Wigner-Seitz model of alcaline metal [11]. The lat- 
ter circumstance turns out to be quite important, since in 
some cases the cavities, where a particle or an atom could 
reside, form a lattice, similar to that of an alcaline metal, 
like certain interstitial sites of a metal supercell, e.g. oc- 
tahedral positions of palladium fee lattice [12]-[14]. In 
this case a particle (or valence atomic electron, provided 
that the whole lattice of cavities is occupied by atoms) 
finds itself in a periodic potential of a cubic lattice, and 
so the description of its ground state could be based on 
the first principles of the Wigner-Seitz model [11]. With 
the same assumptions as in [11], it turns out to be a 
special type of "confinement" under Neumann boundary 
condition in the corresponding Wigner-Seitz cell. 

The purpose of this letter is to explore the features 
of such a type of "confinement" state in a cell, formed 
by a spherical cavity of radius R with an outer potential 
shell of physically reasonable width and depth, and Neu- 
mann condition on the outward boundary. A number of 



nontrivial properties of such state, a part of which being 
similar to those described earlier for atoms trapped endo- 
hedrally inside a fullerene molecule [4] , and more recently 
by means of general reflecting boundaries [15], while an- 
other part being quite new, is discovered by studying 
the asymptotical behavior of energy levels for large R. 
Moreover, such an approach allows for a valuable anal- 
ysis of conditions, under which such phenomena could 
take place. In particular, we describe the case, when the 
ground state of atomic H considered as a function of R, 
contains a deep and strongly pronounced well, where the 
bound energy could be remarkably larger than that of ls- 
level of the free atom E\ s , as well as the situation, when 
the lowest level reveals slowly decreasing power asymp- 
totics for large R and so its bound energy could exceed 
Eis for actual nanocavities with R ~ 100 — 1000 nm. 



2. General treatment of a "not going out" state 

Stationary state of a particle with mass m confined in 
a vacuum cavity Q with boundary £ should be described 
by an energy functional of the following form 



E[ip] 



df 



h 2 -. 



u(f) |v>| : 



(i) 



where U(f) is the potential inside f2, while the surface 
term J" s corresponds to contact interaction of the parti- 
cle with medium, in which the cavity has been formed, 
on the cavity boundary. The properties of this surface 
interaction are given by a real- valued function A(r). 

From the variational principle with normalization con- 
dition (ip\ip) = J n dr \if.<\ 2 = 1 it follows that 
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inside Q combined with boundary condition imposed on 
if) on the surface E 



nV + A(f) V 



, 



(3) 



with n being the outward normal to E. 

Boundary condition (3) is known in mathematical 
physics as Robin's (or third kind) condition, under which 
the spectral problem (2-3) is self-adjoint and so con- 
tains all the required properties for a correct quantum- 
mechanical description of a nonrelativistic particle con- 
fined in fi [15], [16]. The particle "not going out" property 
is fulfilled here via vanishing normal to E component of 
the quantum-mechanical flux 



nj = , 



where 



(4) 



(5) 



At the same time, tangential components of j could be 
remarkably different from zero on E and so the particle 
could be found quite close to the boundary with a marked 
probability. In particular, such a picture takes place in 
the Thomas- Fermi model of many-electron atom [17], as 
well as in quark bag models of hadron physics [18], [19]. 

When A = 0, the interaction of the particle with en- 
vironment is absent and so eq. (3) transforms into Neu- 
mann (second kind) condition 



nVip = , 

E 



(6) 



what corresponds to the boundary condition of confine- 
ment for a scalar field in rclativistic bag models [18]. 
Moreover, condition (6) appears in the Wigner-Seitz 
model of an alcaline metal [11] and describes derealiza- 
tion of valence electrons creating the metallic bond, by 
continuing the atomic wavefunction periodically in the 
lattice. Indeed such a "confinement" state is at the aim 
of our study. 

If A — > oo, then (3) turns into the Dirichlet condition 



V>l s = 



(7) 



and so describes confinement by an impenetrable barrier. 

There are two well-established and quite important in- 
equalities for the ground state energy in the Dirichlet and 
Neumann cases of confinement [16]. The first one takes 
place for the Dirichlet problem (7) and tells, that if the 
volume f2 is embedded in volume f2i, then E(0,) > E(Cti) 
for any nonsingular U (r) . Another one concerns the Neu- 
mann case (6) and gives the following estimate for the 
ground state energy 



EgroundiS^) ^ 



(8) 



The inequality (8) follows immediately from the varia- 
tional principle, if one considers a constant trial wave- 
function in order to fulfill the boundary condition (6) 
in the simplest way. It can be easily generalized to the 
Robin's case (3) in the following fashion. Let us consider 
the confinement state in f2 with no surface interaction, 
but with modified potential function 



U x (f)=U(f) + — A(f)<5 Sl (r), 



(9) 



where S^ir) denotes surface S- function with Ei being a 
surface embedded in CI. The additional term in the mod- 
ified potential U\ gives rise to the following contribution 
to the energy functional 



AE[ip] 



2m 



da A(r) |i/>| 2 



(10) 
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but at the same time doesn't affect the Neumann bound- 
ary condition on E. The latter makes it possible to draw a 
direct analogy with the procedure leading to the inequal- 
ity (8), since the trial wavefunction can be still chosen as 
a constant throughout f2. If we consider now the limit 
Ei — > E from the inside of O, then the region fii C O 
surrounded by the surface Ei tends to Q, while the contri- 
bution of the region — Ox to energy becomes negligibly 
small. Proceeding further this way, we get the expression 
(1) as the limiting point for the energy functional, and 
so the following inequality for the ground state energy 



Egroundi^) ^ 



h 2 J^da A(f) 
2m / n dr 



(11) 



The estimates (8,11) turn out to be quite effective for 
understanding the ground state properties, especially for 
the case of extremely small cavities. 



3. Robin's reflecting boundaries 

Now let us consider the case of Robin's boundary 
condition (3). Since it has been already studied in 
[15], [16], [20], [21], we'll point out here only those details, 
which arc required for dealing with a more complicated 
and realistic model described in the next section. 

First example is a particle in a spherical potential well 
of radius R with a constant potential U(r) = Uq , r < R , 
and surface interaction A =Const [15], [21]. In what fol- 
lows, in order to provide an effective comparison of re- 
sults, obtained for quite different systems, we'll use rela- 
tivistic units H = c — 1, wavenumber and energy will be 
expressed in units of the particle mass m, while distances 
— in units of the particle Compton length 1/m. Consid- 
ering Uq as a reference point for the particle energy, for 
s-levels one obtains 



tan kR 



kR 
1 - XR 



(12) 
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where k = y2E. 

It is easy to see from (12), that the energy levels con- 
sidered as functions of R reveal remarkably different be- 
havior depending on the sign of A. More concretely, when 
A > and so describes reflection between the particle and 
environment, for R — > t he wav enumber of the lowest 
energy level behaves like -\/3A/ R, while the ground state 
energy increases in the following way 



Egroundi.^) ^ on 1 ^ 
ZK 







(13) 



what follows directly from eq. (12) as well as from the es- 
timate (11). Such behavior of E groun d(R) confirms, that 
for confined systems the standard uncertainty relation 
should be replaced by a generalized one, which doesn't 
imply, that for R — > the kinetic energy of the particle 
could be estimated as 0(1/R 2 ) (see [15], [22] and discus- 
sion therein). The latter should be definitely correct in 
the case A — > oo only, i.e. in the case of genuine trapping 
of the particle in a cavity by an impenetrable potential 
wall. For eq.(12) such behavior occurs for the particle 
states with positive energy in the case of surface attrac- 
tion A < 0, when for R —> the wavenumber of the lowest 
positive level behaves like C/R with C = 4.49341 being 
the first root of the equation tanx = x, while the en- 
ergy — like C 2 /2R 2 . For R — s- oo both types of solutions 
for the lowest positive level reveal the same asymptotics 
E{R) — > tt 2 /2R 2 , what corresponds to the Dirichlet con- 
dition (7). 

For A < 0, i.e. for the case of attraction between the 
particle and environment, the generalized uncertainty re- 
lation for confined systems [15], [22] provides, that the 
ground state s-level lyes below the well's bottom and so 
should be found from eq.(12) via k —t in, i.e. from equa- 
tion 



tanh kR 



kR 



1 + \\\R 



(14) 



For R — > the wavenumber n{R) reveals the asymptotics 
y/3\X\/R, hence 



(15) 



what could be easily verified by estimate (11) again. 
There are no contradictions with the general proper- 
ties of the energy spectrum of a nonrelativistic particle 
here, since for A < the surface term in the expression 
(1) could be arbitrarily negative due to \4>\ 2 on the box 
boundary, which might be now arbitrarily large without 
violating the normalization condition. 

It should be noted also, that for R — > oo such a level 
reveals the following asymptotics 



E g round(R) -> -\ 2 /2-\\\/R+0(l/R 2 ) , R 



oo 



(16) 



and so its behavior on the whole half- axis < R < oo 
should be quite similar to a shifted downwards hyperbole. 



Therefore for A < the particle lowest s-level lyes be- 
low the well's bottom even in the case of increasing well's 
radius, but this property cannot be detected from esti- 
mate (11). The latter circumstance should be quite ev- 
ident, since for large R the constant trial wavefunction 
cannot be a good approximation to the genuine wave- 
function of the problem. 

This example shows explicitly, that the spectrum of 
stationary states of a particle confined in a box with 
general "not going out" conditions could reveal features, 
which arc quite different from the deconfinemcnt case. 
In particular, for A < the behavior of the ground state 
is such, that the energetically most favorable state of a 
particle is to be caught by the smallest cavity. 

The "not going out" state of atomic H with nuclei 
charge q in a spherical cavity with radius R and bound- 
ary conditions (3) turns out to be even more specific 
[15], [16], [20], [21]. As in the previous case, surface inter- 
action is given by the constant A, while motionless point- 
like atomic nuclei is in the center of the cavity, then spher- 
ical symmetry is maintained and the ground state en- 
ergy minimized. From the solution of the Schroedinger- 
Coulomb problem for the radial wavefunction of the elec- 
tron state with orbital momentum I one obtains up to a 
numerical factor [17] 



Ri(r) = e" 7 V $(& ; ,c ; ,2 7 r) 



(17) 



where 



7 = V=2E , h = I + 1 - q/<y , ci =21 + 2 , (18) 

and $(&, c, z) is the confluent hypergeometric function of 
the first kind (Kummcr function). Definition, notations 
and main properties of the Kummcr function follow ref. 
[23]. Substituting (17) into the boundary condition (3) 
yields the following equation for energy levels 

[qft + (A - 7 )R - 1] $ R + [I + 1 - qh] = , 

(19) 

where 

<P R = $(&,, Q,2 7 i?), = *(&! + !, cj,2 7 .R). (20) 

As in the previous case of a potential well, the most 
significant changes in the spectrum take place for R — > 0, 
what could be seen at once from the estimate (11). Here 
it should be noted, that for atomic H the limit R — > 
takes some care, since relativistic effects give rise to the 
restriction R > 10 for the cavity sizes, where such an ap- 
proach to the confinement problem, based on boundary 
condition (3), should be valid [21]. So in what follows 
the limit R — > should be understood either as a purely 
mathematical property of equations under consideration, 
or as decreasing R up to R ~ 10. To underline the exis- 
tence of this lower limit, the curves shown on Figs. 2-6 
below will start from R = 10 too. 

There are two types of the lowest level of atomic H 
in dependence on relation between A and q. The first 
one takes place under assumption, that for R — > the 
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wavenumber 7 remains finite, and so in the vicinity of 
R = it could be represented by a series 



E, eV 



7(7?) = 70 + 7i-R + 7 2 i? 2 + 



(21) 



Expanding ,$_r(6+) in a power series in R (what 
is always possible, since the Kummer series converges 
everywhere in the complex plane), to the lowest order 
one obtains from (19) that / = 0, and by proceeding 
further 



\ 22 

A = 1 , 7o = 1 



7„ = 0, n > 1 



(22) 



It follows from (22), that if A = q, then the ground state 
energy of atomic H in a cavity for any < R < 00 
precisely coincides with that of ls-level of the free atom 



E„ 



%d {R) = E ls = -q 2 /2 , 



(23) 



what has been already mentioned in [16], [20]. 

More precisely, for / = 0, A = q, 70 = ±q eq. (19) 
is satisfied for all R. For 70 = q it is provided by 
bo = and $(0, 2,z) = 1, while for 70 = — q one ob- 
tains 6 = 2 , $(2, 2, z) = e z , $(3, 2, z) = (z/2 + l)e z , 
and in both cases substitution into (19) gives an iden- 
tity. There is however no twofold degeneracy of the level, 
since both signs in 70 = ±g correspond to the same ra- 
dial ls-function Ro(r) = Ae~ qr , what should be quite 
obvious, because the parameter 7 is defined via relation 
E = — 7 2 /2, where the sign of 7 isn't fixed. 

As for a particle in a potential well, another type of lev- 
els reveals for R — > asymptotic behavior similar to (13) 
or (15) and is found by assumption, that in the vicinity 
of R = the wavenumber 7 is represented by a series 



7 (J?)= -7=+£o + 6>/fl- 
v R 



(24) 



Substituting (24) into eq.(19), to the lowest order in V-R 
one obtains again I = 0, while higher orders of expansion 
in yield 



q 2 + 3q\ + 6A 2 
20£ 



(25) 

As a result, for such type of s-levels of H in a cavity one 
obtains the following dependence on the cavity radius for 
R -> 

p fmj 3(g-A) g 2 + 3 g A + 6A 2 ^- 

Eground{R) ~> ^ + O(ViJ) , 



2i? 



20 



i? -> 



(26) 



Qualitative explanation of linear dependence on q and A 
is quite simple. As for a particle in a spherical well, for 
R — > the atomic wavefunction of such ls-level inside 
a cavity becomes almost constant, and so the estimate 
(11), which reproduces the first term in (26), turns out 
to be almost exact too. The numerical solution of eq.(19) 



X = q 

X= (1+0.01) 1 
A = (l±0.02)< 



FIG. 1: The lowest s-level of atomic H with q — a in a cavity 
with boundary conditions of "not going out" (3) as a function 
of radius -R. The sign of the shift of A relative to q and 
corresponding shift of levels relative to Ei s coincide. 



for q = a ~ 1/137 and A = (1 ± 0.01)g , A = (1 ±0.02)9 
shows, that the behavior of such s-levels tends to the 
asymptotics (26) for R of order about several tenths of 
a B = l/a~ 137 (Fig.l). 

The analogy between a particle in a well and H in 
a cavity remains valid for R — > 00 too, where it could 
be easily checked by means of asymptotic expansion for 
j ®R.{b+) in (19), that in the case of surface attrac- 
tion A < there exists one more level E{R) with negative 
limiting value E(oo) = — A 2 /2, besides the discrete spec- 
trum of the free atom, and power asymptotical behavior 
for R — > 00 



E(R) 



-\ 2 /2-(q-\)/R + 0{l/R 2 ) , R 



(27) 



For A < — q < this analogy could be extended on the 
whole range of cavity sizes, since under these conditions 
E{R) turns out to be the lowest atomic s-level with the 
form of shifted downwards hyperbole, as for a particle in 
a well. 

Now let us turn to the next type of atomic levels in a 
cavity, which appear under assumption, that jR remains 
finite for R — > 0. To maintain the connection with two 
previous types of levels, we consider only s-levels with 
I = and rewrite (19) in the form 



(2d/dz + X/j-l)^(b,2,z)\ z=2ryR = . 

Since 7 — s- Const/i? for R 0, then A/7 — s- 
1 — q/j —> 1 , and so (28) transforms into 



(2S/a«-l)*(l,2,«)L 



z=2jR 







(28) 
bo = 

(29) 



Taking account of <I>(1, 2, z) = (e z — 1) j z , from (29) one 
obtains 



jR 



tcLIl X 7-; 



(30) 



what describes positive energy levels with the asymp- 
totics 

E n {R)^^, i?-»0, (31) 
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i.e. excited states of a particle (electron) in a well with 
Neumann boundary conditions (6). So all the s-levels be- 
sides Is (provided that the latter turns out to be the low- 
est one and falls down for R — > 0, what implies |A| < q) 
should for R — > reveal asymptotical behavior (31), 
while levels with / ^ lye even higher due to the centrifu- 
gal term. At the same time, for R 1 all the ns-levels 
(as well as levels with I ^ 0) tend to their asymptotical 
values, corresponding to those of the free atom, exponen- 
tially fast 



E n (R) — E„ 



In 



2 A 



A + 7« 



7n#> 1 



where 



~ln/ 2 , In = q/n , 



1,2. 



(32) 



(33) 



are the ns-levels of the free atom. Remark, that lev- 
els with 7„ < A should approach their asymptotics from 
above, while those with 7„ > A from below. 

It should be specially noted, that the asymptotics (32) 
turns out to be an exceptional feature of those confined 
atom levels, which originate from the discrete spectrum 
of the free atom, since such asymptotics is created by 
approaching the argument of the factor T -1 (&), entering 
the asymptotics of the Kummer function $(&, c, z), to the 
pole b — > —n r , n r = 0, 1, ... . Asymptotics for R — >• oo 
of all the other atomic levels in a cavity, which origi- 
nate from the continuous spectrum of the free atom and 
the additional level (27), caused by attractive interac- 
tion with environment, turns out to be a power series in 
1/i?, and their asymptotical values could be either non- 
negative only, or for A < contain one negative point 
E(oo) = -A 2 /2. 

If A = ±7„, the asymptotics (32) modifies in the next 
way. The exponential behavior is preserved, while the 
non-exponential factor undergoes changes in such a way, 
that the ns-levels approach their asymptotics of the free 
atom from above only. For A = 7„ > their asymptotics 
takes the form 



E n (R) - E n -> (n - 1) 



7n 



(2 7 ni?) 



2(n-l) e -2 7 „_R 



7n#> 1 



(34) 



while for the lowest level E\ (R) the exponential part dis- 
appears completely, since in this case A = 71 = q, and 
as it was mentioned above, E\(R) becomes a constant, 
which coincides with E\ s = —q 2 /2. 

For A = — 7„ < instead of (32) one obtains 



E n {R) — E„ 



1 



In 



(2 7 „i?) 2(n+1) e" 



and moreover, the limiting point E(oo) of the level E{R) 
with the power asymptotics (27) coincides with the corre- 
sponding level E n of the free atom (33), what in turn rep- 
resents a remarkable example of von Neumann- Wigner 
avoiding crossing effect, i.e. near levels reflection under 
perturbation [17], [24] — infinitely close to each other for 
R — > 00 levels E n (R) and E(R) should for decreasing R 
diverge in opposite directions from their common limiting 
point E„ . Perturbation in this case is performed by the 
atomic nuclei Coulomb field, since under general bound- 
ary conditions (3) the electronic wavefunction doesn't 
vanish on the cavity boundary, and so for R 1 the 
maximum of electronic density should be shifted into the 
region of large distances between the electron and nuclei, 
where the contribution of the Coulomb field is negligible 
compared to boundary effects. When R decreases, the 
Coulomb field increases, hence E n (R) should go upwards 
according to (35), while E(R) goes downwards according 
to the asymptotics 



E(R) -> E„ 



1 1 
R 



0{l/R 2 ) , R^oo . (36) 



7n-R> 1 



(35) 



So the energy spectrum of atomic H (with q > 0), con- 
fined in a cavity with Robin's condition (3), turns out to 
be the following. For A = q the lowest s-level acquires 
the constant value E\ s of the free atom, for A > — q it 
behaves for R — > according to (26) with an energy shift 
depending on sign (A — q) and for R 1 it approaches 
Eis exponentially fast, while for A < — q < it trans- 
forms into the level E{R) with power asymptotics (27). 
Excited states in all the cases should for R — > reveal the 
behavior (31). And for an H-like atom there once more 
takes place the situation, similar to that for a particle in 
a potential well, namely — whenever A < q, the atomic 
state with largest bound energy, which could sufficiently 
exceed the bound energy of the lowest level of the free 
atom (23), takes place in the smallest cavity. 



4. Atomic H in the Wigner-Seitz cell 

So far, by formulating the confinement problem (2- 
3) it was implied, that a particle in such a "not going 
out" state interacts with environment only on the cavity 
boundary E, i.e. through certain <5-like potential, what 
leads to the surface term in the energy functional (1). 
In a more realistic approach one should consider instead 
of a (5-likc interaction an outer potential shell of nonva- 
nishing thickness d, into which the particle penetrates 
and interacts there with cavity environment. In the limit 
d — > such potential shell should transform into con- 
tact interaction on the surface E. For these purposes the 
boundary condition (3) should be replaced by an equa- 
tion of Schroedinger type, describing particle interaction 
with medium inside the shell, whose potential might be 
quite different from U(r). In the case of spherical cavity 
and shell the first choice for the shell potential is a con- 
stant [/o; as by modelling the cndohedral environment 
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[4]. Then instead of (3) one obtains 
h 2 1 

-— A + U i) = Eif>, R<r < X = R + d , (37) 
2m 

with Neumann condition on the outward shell boundary 
at X = R + d 



dip/dr\ r=x = 



(38) 



There is no A in (38), since the role of interaction with 
environment is played now by eq. (37). Moreover, as 
it was mentioned above, such approach allows for a suf- 
ficiently more wide problem statement, since it doesn't 
imply genuine trapping of a particle by the given volume. 
In particular, boundary condition (38) appears in a quite 
natural way by considering the particle ground state in 
a cubic lattice, formed by cavities of the same type in a 
crystal matrix, within the Wigner-Seitz model [11] with 
the same assumption, that the genuine wavefunction of 
the problem will actually have the highest (0 h ) crystallo- 
graphic symmetry which is not very far from the spherical 
one. The well-known example of such sublattices is given 
by interstitial sites in certain metals and alloys [12]- [14]. 
The cavity together with an outer shell form in this case 
a kind of the Wigner-Seitz cell, while (38) turns out to 
be the condition of periodic continuation of the particle 
wavefunction between neighboring cells. 

It is also worth while noticing, that in the latter case 
instead of the energy level with such periodic wavefunc- 
tion the whole set of states tp^(r) = itg(r) exp (ikf^J with 
u^(r) = u^(f + a), where a is the period of the cavities 
sublattice, while the wavevectors fc fill in the correspond- 
ing first Brilluen zone, should be considered. In any case, 
however, periodic wavefunction describes the level with 
k = 0, hence the position of the bottom of the first Bril- 
luen zone, what is quite important itself. 

Besides this, the magnitude of the shell potential Uq 
should depend on the penetration depth d in such a way, 
that provides a transition to a (5-like interaction with cou- 
pling constant A for d — > 0, what for 



h 2 

U d^ — A 
2 m 







(39) 



Note, that the limits d — > and Uq — > oo don't commute 
— when Uq — > oo with d ^ 0, then one obtains confine- 
ment in a cavity by an impenetrable barrier, hence with 
boundary condition (7), while for d — > and finite prod- 
uct Uod, on the contrary, the case of general boundary 
condition (3) takes place. 

For s-levels of atomic H in such a cavity with an outer 
shell instead of contact interaction, one obtains the fol- 
lowing spectral problem (in units, introduced in sect. 3) 

[kR(1 - KXt&nhKd) + (g/7 - 7R) (tanh nd - kX)] $_r+ 



where k 2 = 2(Uq — E) , while all the other quantities are 
defined as in (18) and (20). It is easy to see, that the 
relation (39) gives 



kR- 



1 — kX tanh nd 
tanh nd — kX 



XR- 1 



, 



(41) 



whence it follows, that for d — > eq.(40) transforms into 
eq.(19) for atomic s-levels with boundary condition (3). 

It should be specially remarked, that the limits d->0 
and R — > don't commute either. In particular, if d ^ 
0, then for the lowest s-level the solution of eq.(40) for 
R — > leads to 



Eground(R) ~> Uq , R -> 0. 



(42) 



The latter could be easily detected from (11), which in 
this case gives 

Eground(R) < E tr ial{R) = 



3R 2 + 3Rd + d 2 
(R + d) 3 



U Q d 



3R 2 



2(R + d) 3 



(43) 



and E tr iai(R — > 0) — > Uq , combined with the above- 
mentioned feature, that for R —> the estimate (11) 
turns out to be exact. 

More precisely, there are two types of solutions of 
eq.(40) for R — > 0. The first one originates from (40) 
by neglecting the term with kR and omitting the com- 
mon factor (tanliKcf— kX) — > (tanhftd— Kd), what leads 
to the following relation 



(q/j - ~fR)$ R + (1 - qft)$ R (b+) = 



(44) 



For R — > eq.(44) contains no solutions with finite en- 
ergy, since when -fR — > 0, then , $_r(6+) — > 1, 
hence (44) reduces to 1 = 0, and otherwise, when 
jR — > Const 7^ 0, then q/j — > and 6o — > 1, thence 
(44) could be simplified up to 



z$(l,2,z) = 2$(2,2,*) 



(45) 



Eq.(45) in turn reduces to e z + 1 = 0, whence 7„ = 
i(ir/2 + nn)/R , what corresponds to a series of highly 
excited s-states with energies 



E„. -> 



(tt/2- 



2R 2 



R^0 . 



(46) 



The second type of solutions of (40) for R — >• emerges 
from the factor (tanned — nd) , what gives n n = ix n /d 
with x n being the solutions of eq. tanx„ = x n , and so 
leads to another series of s-levels, corresponding to the 
energy spectrum of a particle in a well of radius d and 
Neumann boundary condition (6,38) 



+ (1 - g/7)(tanh«;d- nX)>S> R (b+) = , (40) 



(47) 



7 



These levels reveal a finite limit for R — s- 0, while the 
lowest one, corresponding to x = 0, meets the limiting 
value E (R^0) = U . 

It is easy to verify, that there are no solutions of eq.(40) 
for R — > besides (46) and (47). So the effect of infinite 
descent of the lowest level for R — s- 0, which takes place in 
the case of contact surface interaction for A < q, doesn't 
occur for the potential shell of nonvanishing width. 

The physical meaning of series (46,47) should be quite 
clear. The levels (46) correspond to s-states of continuous 
spectrum of the free atom, when the latter is confined in 
a cavity with R — > 0, while the levels (47) originate from 
ns-levels with exponential asymptotics (32) and a finite 
number of levels E k with power asymptotics for R —> 
oo, which appear for Uq < and turn out to be direct 
analogies of E(R) for the case of contact interaction with 
A < (27) . 

Compared to the case of contact interaction (32), the 
asymptotics of ns-levels for R 1 is modified in the 
following way 

E n (R) — E n — > 



In 
111 



tanh (n n d) - 7.,, 

Kin 

tanh (n n d) + 7,, 



7n#> 1 



{2 ln R) 2n e~ 2 ^ R 



where 



n n = yj2U + 7 



2 

n ' 



(48) 



(49) 



while E n and j n are defined as in (33). It follows from 
(48), that for 



|k„ tanh (k„g!) \ < 7,, 



(50) 



the curves E n (R) approach the ns-levels of the free atom 
(33) for R ^> 1 from below, while for |k„ tanh (K n d) | > 
7„ from above. Therefore the curves E n (R) could for 
finite R reveal nontrivial minima, which lye below the 
corresponding ns-levels of the free atom (33), provided 
that the relation (50) is satisfied. The specific feature of 
the problem with an outer shell is that now such a mini- 
mum, and the deepest one, exists for the lowest s-level as 
well, whereas for vanishing width of the shell and A < q 
this level should for R — s- reveal an infinite falldown, 
and so nontrivial minima could appear for excited states 
only. A crude estimate for such a minimum for the low- 
est s-state can be received from inequality (42) by solving 
dE tna i(R)/dR\ Ro = 0, what gives 



Rn — 



2qd 



2U d 



(51) 



and 



1 2U d - q 
d (2U d-3qf 



4 (U d) 3 - 16q {U df + 19q 2 U d - 6q 3 



(52) 



The cq. (51) predicts the existence of a nontrivial min- 
imum for the lowest state only when 2U^d < q, what is 
more crude, than the exact relation (50). The difference, 
however, should be quite clear, since the estimate (11) 
works well only for small cavities of such type, hence for 
small d, when 2Uod should be identified with A and so 
2Uod < q is nothing else, but the relation A < q. 

As for the boundary condition (3), the asymptotics 
E k (po) = Ek of power levels E k {R) with negative limiting 
values for R — > 00 is found from (40) by taking account 
of the main exponential term in the asymptotics of the 
Kummcr function, what yields the following relation 



k k tanh (k k d) = 



(53) 



where k k = ^/2U Q + % , E k = -7^/2. Note, that if 
Ik = 7n, i-e. the levels E k and E n possess the same 
limiting value for R — > 00, then the l.h.s. of (53) coincides 
with the denominator in the asymptotics of exponential 
levels (48). So vanishing denominator in (48) implies 
once more the change in the asymptotical behavior of the 
exponential level due to the Neumann- Wigner reflection 
effect, what is discussed in detail for the case of the lowest 
level below. 

It follows from (53), that such power levels with E k < 
might appear only for Uq < 0, when (53) takes the form 



7 2 tan 



(v/2|t/ |-7 2 d) = 7 , (54) 



and is nothing else but the equation for even levels in 
one-dimensional square well of width 2c? and depth Uq. 
Therefore the levels E k exist for any Uo < and d > 0, 
their values lye in the interval Uo < E k < 0, while their 
total number K is defined from ir(K — 1) < 2|Z7o|g? < 
ttK , K=l,2,.... 

The asymptotics of the levels E k (R) for R — > 00 takes 
the form 



q 



E k (R)^E k -^- 



\U a 



E k 



1 



Ik 



R \U \(l + %d) R l + %d 



R 



-0{l/R 2 ) 



(55) 



Note, that in such a problem there exist other power 
levels with E k > , which correspond to imaginary 
wavenumbers 7^. and so should be found from the asymp- 
totics of the Kummer function including the power term 
besides the exponential one, but these levels lye wittingly 
higher, than the power (55) and exponential (48) ones, 
whereas our main interest is first of all bent on the lowest 
atomic levels. 

Another crucial difference between power E k {R) and 
exponential E n (R) levels is that the origin of the formers 
is the attractive interaction between the particle (atomic 
electron) with cavity boundary (outer shell) , rather than 
the interaction with the inner shell (atomic nuclei). In 
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FIG. 2: The lowest ls-level of atomic H with q = a as a 
function of cavity radius R for Uo = 10 eV and d = xas- 



FIG. 3: The lowest ls-level of atomic H with q = a as a 
function of cavity radius R for Uo = — 10 eV and d = xcib- 



fact, (53) is equivalent to the equation, denning the s- 
levels of a particle in the attractive potential of a spher- 
ical shell 

U(r) = U 6{R<r < R + d) (56) 

with Neumann boundary condition (38) on the outward 
boundary, in the limit R — > oo. In a slightly different 
language, this circumstance has been pointed out in [4, 
15]. 

A more detailed analysis of eq.(40) turns out to be 
most conveniently performed by means of its numerical 
solution for a concrete set of parameters Uq and d, cor- 
responding to realistic scales of microcavities, in which 
such a "confined" H state could occur (to simplify the 
discussion, henceforth we'll deal with energy in eV). For 
\Uq\ this is 1 — 100 eV, for d — fractions of the Bohr 
radius = 1/a ~ 137, more concretely d — xas with 
x = 2~ p , p = — 1, 0, 1, 4, where the largest 2cib is 
chosen according to the mean width of one-atom surface 
shell, while the smallest as/16 — in accordance with the 
lower limit, following from rclativistic effects [21]. The 
range of values for \Uq\ is defined by taking into account, 
that | Uq I could vary from ~ 1 eV for vacuum "bubbles" in 
supcrfluid He 4 [6] up to dozens eV in quantum chemistry 
[l]-[5]- 

The most simple and transparent example is given 
by the potential barrier Uq > 0. In this case the 
lowest energy atomic H state coincides with the expo- 
nential ls-level, whose behavior as a function of R for 
Uq = 10 cV is shown on Fig. 2 (for R > 10). In ac- 
cordance with relation (50), which in this case gives 
d < 100, there are pronounced minima for such Uq and 
d = a^/4 , as/8 , as/16 on the curves E±(R), and so a 
cell with such parameters turns out to be an effective H- 
trap. A minimum exists also for d = as/2, but here it is 
quite weak (bound energy exceeds only ~ 13, 9 eV) and 
takes place at R = 289, while for d = clb , 2as, on the 
contrary, the H state with the lowest energy is achieved 
for infinitely increasing R. 

The behavior of the lowest level reveals a more pro- 



t/ocrit, eV 



1(1!) 200 300 400 500 




FIG. 4: The "critical" potential U cr ; t as a function of d. 

nounced dependence on d, as well as on Uq, in the case 
of attraction in the outer shell (Uq < 0) due to increased 
amplitude of electronic wavefunction near the outward 
boundary. In particular, for Uo = — 10 cV (Fig. 3) the re- 
lation (50) is fulfilled for all d, hence nontrivial minima in 
the bound energy exist now for all d, including d = 2a_B, 
when a minimum with bound energy 14, 8 eV is achieved 
at R = 254. 

Now let us consider the most interesting case of energy 
levels reconstruction for atomic H in a cavity with such 
an attraction in the outer shell, that the lowest atomic 
level turns out to be the power one Ei(R). For these 
purposes the "critical" potential Uq cr it, which provides 
the coincidence of £4(00) with the limiting value of the 
first exponential level £4(00) = E\ s , what implycs 

7i = 7i , (57) 

should be firstly determined. 

The form of U cr it as a function of the shell width d 
is shown on Fig. 4. Above Uq cr it(d) there lyes the region 
of Uq and d, where the lowest level turns out to be the 
exponential E\(R), while below — the region, where the 
power Ei(R) is the lowest. 
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It is easy to find from eqs.(53, 57), that for d — > oo 
the limiting value of U cr it(d) should coincide with the 
lowest level of the free H (23) with the following asymp- 
totics 



E, eV 



U, 



cnt 



8d 2 4gd 3 



+ 0(l/d 4 ), d^oo 



(58) 

while for d — > the "critical" potential decreases accord- 
ing to 



U crit(rf) -> -q/2d , d -> . 



(59) 



The numerical values of C/ C rit f° r d under consideration 
are presented in Tab.l (besides d = as/16, since in this 
case \U cr it | turns out to be too large). 



x = d/as 


1/8 


1/4 


1/2 


1 


2 


4 


16 


U CTit' eV 


-118.9 


-64.2 


-37.0 


-23.7 


-17.6 


-15.0 


-13.7 



TAB. 1: The values of U 

crit f° r d — xas- 

The behavior of the curves Ei(R) and E\{R) for {/ = 
0.9 t/ cr it and Uq = Uq CT n is shown on Fig. 5 for d = as- 
It is easy to recognize here the effect of avoided crossing, 
discussed in [4, 15], which shows up now in the change 
of asymptotics of the lowest level, when the value of Uq 
coincides with Uq cr ;t. 

It should be emphasized, that the change of the lowest 
level asymptotics for R — ¥ oo from exponential into the 
power one takes place indeed when Uq reaches Uq ct h 
from above, not earlier and not later. In this case 
the limiting point E\ s of the free atom is the same for 
the curves E 1 (R) and Ei(R), and so as for the case of 
contact interaction (35,36), the exponential level Ex(R) 
should approach its limiting point from above, thence the 
power Ei(R) turns out to be the lowest one due to the 
Neumann- Wigncr reflection. Let us underline specially, 
that it is indeed an exchange of asymptotical behaviour 
for the lowest level — the levels E\ (R) and E\ (R) could 
be infinitely close to each other, but don't touch and 
all the more that intersect, since all the s-levels in such 
a problem cannot be degenerate. Note also, that this 
exchange of asymptotical behavior proceeds in the fol- 
lowing way — when Uq Uq cr it from above, the de- 
nominator in the r.h.s. of (48) tends to zero, hence the 
exponential tail of the curve E\ (R) is shifted to more and 
more large R, for Uq = Uq cr jt it disappears completely, 
and so the power behavior extends to the whole half-axis 
< R < oo. 

The behavior of the level Ei(R) for Uq = Uq C rit(^)) i- c - 
at the moment when it becomes the lowest one, is shown 
on Fig. 6 for d = 2~ p as with p = —1,0, 3. (The case 
d = as/16 is omitted, since Uq and the upper limit of 
bound energy acquire in this case too large values Uo ~ 
—240 cV, bound energy > 300 cV, and so in background 
of the curve for d = as/16 the details of the other curves 
become illegible.) 



U = 0.9 C/ocit 



Ei(K) 
Et(R) 



E, eV 



Uo = £/ 0c 



Ei(R) 



FIG. 5: Ei(R) ,E!(R) for U = 0.9 U crit and U = U crit 
for d = as- 



£,(R),eV 



200 300 400 500 



x = 2 
x = 1 
■ x = 1 / 2 
x= 1/4 
x = 1 / 8 



FIG. 6: The behaviour of E^R) for U = U crit (d) and 
d = iob. 



Fig. 6 shows explicitly the power asymptotics of the 
level Ei(R) for R — > oo. The shift of power levels rel- 
ative to their asymptotical value decreases sufficiently 
slowlier than that of the exponential ones, which arrive 
at their asymptotics for R of order of several as already, 
and so such a "confined" atomic H state turns out to 
be energetically favorable up to actual nanoscales, pro- 
vided that Uq < Uq cr it(d). Numerical values (in eV) 
for the shift of bound energy of the power level E\(R) 
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relative to the free H are given in Tab. 2 for a cavity of 
nanosize with Uq = U C vit(d), i-e. at the moment when 
Ei(R) transmutes into the atomic H ground state, for 
d = a B , as/2, a B /4. 



R 


1 nm 


10 nm 


100 nm 


1000 nm 


AE{d = 


= a B ) 


1.00894 


0.102517 


0.0102667 


0.00102686 


AE(d = 


as/2) 


1.5722 


0.157199 


0.0157202 


0.00157203 


AE(d = 


as/4) 


2.08505 


0.20661 


0.0206425 


0.00206407 



TAB. 2: The values of AE = E Xs - E 1 (R) in a nanocavity 
with Uo — U cr n(d), when Ei(R) transmutes into the lowest 
atomic level. 



5. Conclusion 

To conclude let us firstly mention, that a single cavity 
might be just a simple hollow cage without any special 
confining property, but a large set of cavities, forming a 
cubic lattice, could reveal such properties due to quantum 
coherence effects, similar to those creating the metallic 
bond in the Wigner-Seitz model [11]. As a result, each 
single cavity of such a lattice transforms into a kind of 
the Wigner-Seitz cell, formed by a cavity with an outer 
potential shell. 

The properties of the particle state in such a cell turn 
out to be quite different from confinement by a potential 
barrier [9], [10]. An example of such kind is presented 



in sect. 4 by atomic H in a cavity with potential shell, 
provided that the set of cavities, occupied by atoms, 
forms a crystal structure similar to that of an alcaline 
metal. In particular, in dependence on the outer shell 
parameters the upper limit for the bound energy of H in 
such a cell could be more large, than several times the 
bound energy of the lowest ls-level of the free atom. At 
the same time, in the case of a power lowest level the 
bound energy decreases very slowly for increasing cavity 
size, therefore such a state should be energetically favor- 
able compared to the free atom up to actual nanocavi- 
ties with R ~ 100 — 1000 nm. The latter circumstance 
means, that artificial macroscopic lattices, created from 
such nanocavities in suitable media, could serve as quite 
effective containers of H. Even more interesting results 
for searching possible new effects appear in the case of 
more complicated atoms and simplest diatomic molecules 
in such a cell and a lattice with the same parameters as 
employed in sect. 4 [25]. 
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